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Correction to the hyperfme constant A of an external s-electron in many-electron atom caused 
by the Breit interaction is calculated analytically: SA/A = 0.68Za 2 . Physical mechanism for this 
. correction is polarization of the internal electronic shells (mainly Is 2 shell) by the magnetic field 

' of the external electron. This mechanism is similar to the polarization of vacuum considered by 

Karplus and Klein Q long time ago. The similarity is the reason why in both cases (Dirac sea 
polarization and internal atomic shells polarization) the corrections have the same dependence on 
' the nuclear charge and fine structure constant. 

In conclusion we also discuss Za 2 corrections to the parity violation effects in atoms. 

(N . 
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I. INTRODUCTION 



Atomic hyperfme structure is caused by the magnetic interaction of the unpaired external electrons with the nuclear 
magnetic moment. There are two types of relativistic corrections to this effect. The first type is a single particle 
correction caused by the relativistic effects in the wave equation of the external electron Q . This correction is of the 
order of (Za) 2 , where Z is the nuclear charge and a is the fine structure constant. Solving the Dirac equation one can 
find this correction analytically in all orders in (Za) 2 [||. An alternative way to find this correction is direct numerical 
solution of the Dirac equation. Correction of the second type has a many-body origin: it is due to polarization of 
paired electrons by the Breit interaction of the external electron. There are two kinds of paired electrons in the 
problem: a) Dirac sea, b) closed atomic shells. The contribution related to the Dirac sea was calculated by Karplus 
and Klein long time ago [|): SA/A = a/2n — Za 2 (5/2 — In 2). The a/2ir part is due to usual anomalous magnetic 
moment and Za 2 part comes from a combined effect of the Nuclear Coulomb field and the Breit interaction. Wc 
would like to stress that as soon as Z > 12 the Za 2 part is bigger than a/27r . Effect of the polarization of atomic 
shells by the Breit interaction has been recently calculated numerically 0-0. These calculations were performed for 
Cs (Z = 55) because they were motivated by the interest to parity nonconservation in this atom. Results of these 
calculations are somewhat conflicting, but nevertheless they indicate that the correction for an external s-electron is of 
. the order of ~ ±0.4%. In spite of being rather small this correction is comparable with the present accuracy of atomic 
many-body calculations and therefore it must be taken into account. The purpose of the present work is to calculate 
analytically the correction induced by the Breit interaction. This allows to elucidate the physical origin of the effect 
and its dependence on the atomic parameters. This also provides an important lesson for a similar correction to the 
O , parity non-conserving amplitude which we discuss in the conclusion. 

>>, 

^ . II. CONTRIBUTION OF THE DIRECT ELECTRON-ELECTRON MAGNETIC INTERACTION 

Oh; 

In the present work we do not consider single particle relativistic effects (Dirac equation), so we assume that in 
zero approximation the atom is described by the Schroedinger equation with Coulomb electron-nucleus and electron- 
electron interaction. For magnetic interaction throughout this work we use the Coulomb gauge. Vector potential and 



X 

magnetic field of the nucleus are of the form 



Aiv(r) = (1) 



87T 

HN(r) = V x A w -> —^ N S(r), 

where /iN is the magnetic moment of the nucleus. We keep only the spherically symmetric part of the magnetic field 
Tin because in this work we consider only s-electrons. Interaction of the magnetic moment [i\ of the external electron 
with nuclear magnetic field gives the hyperfme structure 
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H hfs = (-mi • H N ) = -C (pn ■ ntf) = A (s • I), (2) 



A = ^^C. 



s I 

Here i/j e (f) is the wave function of the external electron, and A is the hyperfine constant. 
Vector potential of the external s-electron is of the form 

A e (r) = J Ml ^ r ^ 3 R) iPUR)d 3 R = -Mi x V r J |^ d 3 i? = 4tt ^ e 2 (#) i? 2 ^- 

Hence the magnetic field is 



(3) 



W e (r)=VxA e -,yi 2 (r) ft . (4) 

Let us repeat once more that we keep only the spherically symmetric part of the magnetic field. 

The Hamiltonian of an internal electron in the magnetic field A = A^v + A e , Ti = TLn + W e is given by 

(v - -A) 2 

g = c - M2 ■ H + U(r), (5) 
2m 

where ^2 is the magnetic moment of the internal electron. Certainly |mi| = = 27^' nowever directions of fix 
and ^2 are independent. Having eq. (||) in mind one can easily draw diagrams describing correction to the hyperfine 
structure due to the electron-electron magnetic interaction. This diagrams are shown in Fig.l 
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FIG. 1. Diagrams describing the direct Breit interaction correction to the hyperfine structure. The diagrams a) and b) give 
the paramagnetic contribution, and the diagram c) gives the diamagnetic contribution. N denotes the nucleus, e denotes the 
external electron, n denotes the internal electron, and finally m denotes a virtual excitation of the internal electron. The wavy 
line shows magnetic interaction with the external electron, and the dashed line shows magnetic interaction with the nucleus. 

Two equal paramagnetic contributions are given by the diagrams shown in Fig. la, b. Corresponding energy correc- 
tions are 

SE a = 5E b = H(-M2-tt e |#„>, ( 6 ) 

n£ filled 



where 

c, \ - H(-M2 - n N )\n) 

Sip n = 2^ B B H- ( 7 ) 

The diamagnetic contribution shown in Fig.lc is given by the A 2 term from the Hamiltonian (||), hence 

e 2 

5E C = -{n\A e - A N \n). (8) 

mr 

Before proceeding to the accurate calculation of 8E it is instructive to estimate a magnitude of the correction. Let 
us look for example at the diamagnetic correction (g). According to eqs. (0), (^) the vector potentials are An ~ [in jr 2 
and A e ~ fiirip 2 (0). Hence the correction (w) is of the order of SE ~ <? hi^in4> 2 (0)/ (mc 2 r), where r is the radius 
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of the internal shell. Since all the interactions are singular it is clear that the main contribution comes from the 
K-shell, so r ^ clb / Z (as is the Bohr radius). Together with eq. (|j) this gives the following relative value of the Breit 
correction to the hyperfine constant 5C/C ~ Za 2 . So we see that this effect has exactly the same dependence on the 
atomic parameters as the Dirac sea polarization considered in the paper jlj . 

Now let us calculate the coefficient in the Za 2 correction. We consider explicitly only Is and 2s closed shells 
and we also need to consider the external s-electron. In atomic units the single particle energies of these states are: 
Ei = —Z 2 /2, E2 = —Z 2 /8, E e w 0. At small distances the nuclear Coulomb field is practically unscreened and hence 
the wave functions are of the simple form S 
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1p2 



^=e-"/ 2 (l-p/2), (9) 

/87T 



Here p = Zr and J\{x) is the Bessel function. The functions ^1,2 are normalized in the usual way: J ip 2 d 3 p = 1. 
The wave function of the external electron is normalized by the condition if> 2 (0) = 3/(87r). With this normalization 
the leading order hyperfme constant (^|) is equal to unity, C = 1, and therefore this normalization is convenient for 
calculation of the relative value of the Breit correction to the hyperfine constant. Using eqs. (0) , (f|) , (||) an d performing 
summation over spins in the closed shells one finds the diamagnetic correction 



SE C = Za 2 {pi ■ p N ) y 



(^^ d3p / P ^ (pW ) =°- 230Za2 ^-VN). (10) 



The numerical coefficient was found by straightforward numerical integration. Contributions of the inner shells drop 
down approximately as 1/n 3 , so 0.230=0.207+0.023, where the first contribution comes from the ls-shell and the 
second contribution comes from the 2s-shell. 

To calculate the paramagnetic contributions (||) we use corrections Sipn defined by eq. (^) and calculated in the 
Appendix. Substitution of (^), (Q), and ( |A9| ) into formula (|6|) and summation over electron spins in the closed shell 
gives the following result 

8 / f°° 1 f°° \ 

5E a +5E h = --Za 2 (p 1 -p N )I^J eT 2 ' J 2 (y/^) Wl (p)dp + - j e -"(l - p/2)J 2 (V%P)Mp)dp) (11) 

g 

= --Za 2 (/ii ■ ^jv)(0.219 + 0.021) = -0M0Za 2 (pi ■ p N ), 
o 

where we present explicitly the contributions of Is- and 2s-shells. The numerical coefficient is found by numerical 
integration. Similar to the diamagnetic term the contributions of the inner shells drop down approximately as 1/n 3 . 

The leading order hyperfine structure is given by eq. (^|) with constant C = 1 due to the accepted normalization. 
According to eqs. (|l^) and (|ll|) the total correction caused by the direct magnetic interaction is SE a + SEi, + 6E C = 
—QA10Za 2 (pi ■ /1 at). Comparing this with eq. @) one finds the relative value of the direct correction: 

6 -± r - = 6 \- = OAlOZa 2 . (12) 

III. CONTRIBUTION OF THE EXCHANGE ELECTRON-ELECTRON MAGNETIC INTERACTION 

Exchange diagrams contributing to the correction are shown in Fig2. 
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FIG. 2. Diagrams describing the exchange Breit interaction correction to the hyperfine structure. The diagrams a) and b) 
give the paramagnetic contribution, and the diagram c) gives the diamagnetic contribution. N denotes the nucleus, e denotes the 
external electron, n denotes the internal electron, and finally m denotes a virtual excitation of the internal electron. The wavy 
line shows magnetic interaction with the external electron, and the dashed line shows magnetic interaction with the nucleus. 

The diagrams Fig. 2a, b show the "paramagnetic" contributions, and the diagram Fig. 2c shows the "diamagnetic" 
contribution. Note that the contributions of the diagrams Fig. 2a, b must be doubled because the opposite order of 
operators is also possible. Let us begin with the "dimagnetic" term. Comparison of Fig.lc and Fig. 2c shows that the 
direct and the exchange contributions are very similar and therefore the simplest way to derive the exchange term 
is just to make appropriate alterations in eq. ([lO]) which gives the direct contribution. The alterations are obvious: 
l)opposite sign, 2)i/>„ — > ip e i/) n , 3)V> 2 ~~ * V'eV'nj 4)there is no summation over the intermediate spins, hence 4/3 — ► 2/3. 
Thus the result is 

SEW = -Za 2 (p, (§ - iP) p t d3 P [ ^ We (j/)dV ) = -0.107ZoV ' m)- (13) 

The coefficient is found by numerical integration: 0.107=0.096+0.011, where the first contribution comes from the 
ls-shell and the second contribution comes from the 2s-shell. 

The paramagnetic exchange contribution shown in Fig. 2b is similar to the direct ones given by Fig. la, b. The only 
difference is in algebra of Pauli matrixes and in additional sign (-). This consideration shows that the paramagnetic 
exchange contribution is equal to half of that given by eq. ( |ll| ) 

A / />OC -1 />00 \ 

§ E^ = --Za 2 (/ii • Pn) U e- 2 "J!(y/Tp)w 1 (p)dp + - J e-"(l - p/2)j!(^/8~p)w 2 (p)dpj (14) 
= -0.320Za 2 (>i • Pn)- 



Note that the sign of the exchange contribution is the same as the sign of the direct one (|11|) . 

The diagram shown in Fig. 2a does not have analogous direct diagram because it has the hyperfine interaction 
attached to the line of the external electron. Nevertheless the calculation of this diagram is quite similar to the 
calculation described by eqs. (Q) and After substitution of Syj e from ( |A13| ) and performing summation over the 
polarizations inside the closed shell one finds the following expression for the diagram shown in Fig. 2a 



/ roc 1 f° 

= SirZa 2 ^ ■ p N ) U e' 2p J 1 (y/8^)N 1 (y/8^)pdp+-J 
= 0.156Za 2 (fii ■ pn)- 



p(1 - p/2) 2 J 1 (^)Ar 1 (^)prfp ) (15) 



The total exchange magnetic correction is SEa^ + SE^' + 8E[ ex ^ = —0.271Za 2 (pi ■ /ijv). Comparing this with 



b . 

eq. (H) one finds the relative value of the exchange correction: 



5A^ SC^ 9 , s 

— ^— = = 0.271Za 2 . (16) 



IV. TOTAL BREIT CORRECTION. 
Za 2 CORRECTION DUE TO ELECTRON-ELECTRON COULOMB INTERACTION 

Adding the direct (|l^) and the exchange (|l6|) contributions one finds the total Breit correction 

M£ = ^£ =0 .681Za 2 . (17) 

In the calculation we have not used the explicit form of the Breit interaction, but nevertheless this is the correction 
generated by the interaction which reads in the relativistic form and in the Coulomb gauge as (see ref. ) 

tt 1 ( (ai ■ r)(q 2 ■ r) ^ 

H B = -— ■ a 2 + - 2 j. (18) 

Here r = ri — r 2 is distance between the electrons, and is the a-matrix of the corresponding electron. The Breit 
interaction correction to the hyperfine structure of Cs was previously calculated numerically in the works 



4 



but results of these calculations were somewhat conflicting. Our result ( |l7| ) agrees with that of the most recent 
computation @. Note that eq. (|l7| ) gives the leading in Z part of the Breit correction. There are other parts, say the 
correction to the energy of the external electron which directly influence the hyperfine constant. However the other 
parts contain lower powers of Z . 

The correction (|l^) does not include all Za 2 terms. To realize what is left let us look at the electron-electron 
interaction Hamiltonian in (v/c) 2 approximation || 



H B =a 2 { - 



7r<5(r) - — (Sep + r -^f) pi a P2!3 + ^3 (-(si + 2s 2 ) • [r x Pl ] + (s 2 + 2si)[r x p 2 ]) (19) 
si • s 2 3(si • r)(s 2 • r) 8?r N 

+ ~ ~5 ' 82 ^ 

Here Pi,Si denote the momentum and the spin of the electron. All the terms containing momenta vanish for s- 
electrons, the two last terms are already taken into account by the calculation performed above, however the first 
term has not been considered yet. The matter is that in spite of being a (u/c) 2 -correction it has a nonmagnetic origin. 
It comes from the (v/c)-expansion of the electron-electron Coulomb interaction -(u\ui)(u 2 U2), where it.; is the Dirac 
spinor of the corresponding electron. This is why this term is accounted automatically in the Dirac-Hartree-Fock 
calculations |ic|,^-|7|] . Nevertheless if one wants to separate the total Za 2 correction analytically then the first term 
must be also considered explicitly. Since this term is spin independent, it can contribute only via exchange diagrams. A 
straightforward calculation very similar to that performed above gives the following result for the Coulomb correction. 

SEcouiomb = 2Za 2 {p 1 ■ p N ) (J e- 2p J 2 (y/8^)w 1 {p)dp + 1 J e -"(l - p/2)J*{^8p~)w 2 (p)dp^ 

-AwZa 2 ^!-^) Qf e- 2 <>J 1 (y/8p')N 1 (y/8p')pdp+^ J e' p {l - p/2) 2 J 1 ( y /8p')Ni(y/8p')pdp^ (20) 
= 0.558Za 2 (/ii • pn)- 
Comparing this with eq. @ one finds the relative value of the Coulomb correction: 

= — — — O.bbsZa . (21) 

AC 



Combination of the magnetic (17) and of the Coulomb (21) corrections give the total Za 2 correction to the hyperfine 



constant due to the polarization of the closed electronic shells 

SA 



^ =0.123Za z . (22) 



V. BREIT INTERACTION CORRECTION TO THE PARITY NONCONSERVATION EFFECT 

The Breit interaction correction exists for both nuclear spin independent (weak charge) and nuclear spin dependent 
(anapolc moment) weak interactions. The correction to the spin independent effect is the most interesting one because 



of the high precision in both atomic theory (T^| and experiment 12 ljj. Recent computations ||[7| show that the 
atomic Breit correction to the nuclear spin-independent parity nonconservation (PNC) effect in Cs is about 0.6%. 
This is enough to influence interpretation of the experimental data, see Refs. [|l3]||. The Breit correction to the PNC 
effect can be calculated analytically similar to the hyperfine structure correction. However this calculation is out of 
the scope of the present work. In the present paper I just want a) to estimate this correction parametrically, b) to 
comment on the importance of the Dirac sea polarization. Let us first estimate the correction. It is given by the 
diagrams similar to that shown in Fig. 2a, b. The only difference is that the electron-nucleus magnetic interaction must 
be replaces by the electron-nucleus weak interaction |l3| 

H W = — %-Qw[s-pS(r)+5(v) S -p], (23) 
2y 2m 

where G is the Fermi constant, and Qw is the weak charge. The relative value of the Breit correction is 
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2/3 2 ? Z 

S(H W )/(H W )~^ zr~Za\ (24) 

Here a 2 jr is the magnetic interaction between external and internal electrons, and AE is the excitation energy of the 
virtual state of the internal electron. For the estimation we take K-electrons, therefore r ~ 1/Z, and AE ~ Z 2 . Thus 
the Breit correction to the PNC effect has exactly the same dependence on atomic parameters as the Breit correction to 
the hyperfme structure. It is known that there is a large relativistic enhancement factor for PNC effect p^ |, therefore 
one might think that the nonrelativistic expansion used in the estimate ( pi|) has very poor accuracy. However it is 
not the case, the matter is that the large relativistic factor appears from the distances r ~ nuclear size -C ag/Z, 
therefore this factor is more or less the same for the external electron and for the K-electron. So it is canceled out in 
the ratio ( pi] ) . By the way a similar argument explains a relatively high accuracy of the nonrelativistic expansion for 
the hyperfme structure correction (p"7|). 

I would like also to comment on the importance of the Dirac sea polarization. The effect calculated in refs. HQ and 
estimated in eq. (|24|) is caused by the polarization of the closed atomic shells. This is analogous to the contribution 
( |l7| ) for the hyperfme structure. However one has to remember that for the hyperfme correction there is also the 
contribution of the Dirac sea polarization [El, a/2n — Za 2 (5/2 — In 2) a/2n — 1.81Za 2 , which is bigger and has 
the opposite sign. Only account of both effects together (atomic shells + Dirac sea) has the physical meaning. The 
radiative correction to the nuclear weak charge is known in the single loop approximation, see Ref. jT5| . This includes 
a/ir and a/irln(Mz/n) terms (Mz is the Z-boson mass and /i is some infrared cutoff). However Za 2 radiative 
correction to the weak charge has not been considered yet and it is quite possible that this correction is larger than 
the a /ir contribution (at least this is so for the hyperfme constant). To calculate Za 2 radiative correction one needs 
to go to the two loop approximation, or to work in the single loop approximation but with the Green's functions 
in the external nuclear Coulomb field. Thus: 1) account of the Breit correction calculated numerically in Refs. ||j7]] 
and estimated in eq. ((5J) without account of the Dirac sea contribution is not sufficient, 2) account of the Dirac sea 
polarization can influence agreement between theory and experiment. 



VI. CONCLUSION 



Correction to the hyperfme constant A of an external s-electron in many-electron atom caused by the polarization 
of inner atomic shells by the electron-electron Breit interaction is calculated analytically: SA/A = 0.68Za 2 . This 
correction has the same origin as the Dirac sea polarization effect 5 A/ A — a/2ir — Za 2 (5/2 — In 2) calculated by 
Karplus and Klein long time ago . 

It has been shown that the parametric estimate for the Breit correction to the parity nonconservation effects is also 
Za 2 . We stress that to take this correction into account one needs to consider both polarization of the inner atomic 
shells and polarization of the Dirac sea. 
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APPENDIX A: CORRECTIONS TO THE ELECTRONIC WAVE FUNCTIONS DUE THE HYPERFINE 

INTERACTION 



Correction Sip to the single particle wave function is given by eq. (|7|). Using this formula one can easily prove that 
Sip satisfies the following equation 

{H - E n )Si> n =-J2 \m)(m\V\n) = -V\n) + (n\V\n)\n). (Al) 

Here 

V = -fi-H N = -^-(fx- t i N )S(r) (A2) 

is the perturbation and 



G 



is the Hamiltonian of the Coulomb prob lem. The screening of the Coulomb field is neglected because we consider only 
small distances, r ~ 1/Z. The eq. (Al) has an infinite set of solutions. To find the correct one we have to remember 
that there is the additional condition of orthogonality 



{Sip n \tp v 



0, 



which follows from cq. (|7|). Having the perturbation (A2) it is convenient to represent 6ip n in the form 



(A4) 



(A5) 



where p = Zr. Substitution of ( |A5| ) into ( Al) shows that the functions f(p) obey the following equations 

2 



Is 

2s 



A p + --l)f 1 (p) = -4w5(p) + 4e-P, 



h(p) 



-^8{p) + \e-fl\\-pl2). 



2 1' 

To satisfy the boundary conditions at p = and atp = oo it is convenient to use another substitution 



/i = 



1 



Vl(p), 



(A6) 



(A7) 



h = -e 
P 



-p/2 



where Wi(0) = 1 and Wi(p) grows at larg e p not faster than a polynomial. Straightforward solution of eqs. (A6) 
together with the orthogonality condition ( ]A4[ ) gives 



w 1 = 1 - 2p[ln(2p) - 5/2 + c] - 2p 2 , 

w 2 = 1 - 2p[(l - p/2) Inp- 3/4 + c] - (13/4 - c)2p 2 + p 3 /A, 



(A8) 



where c = 0.577215 is the Euler constant. Altogether eqs. (|A5|) , (|A7|), and dA8|) give 

5^ = J-Z 5 l 2 {p-p N )-e-» Wl {p), 
3a/7t p 

S^ 2 = -^—Z^l\p ■ p N )-e-^ 2 w 2 (p). 

3V87T p 



(A9) 



We also ne ed t o know the hyperfine correction to the wave function of the external electron. Basically it is also 
given by eq. (Al), but t here is a special point concerning normalization. We use an artificial normalization condition 
tpl (0) — 3/ (8tt), see eq. (|A8|). On the other hand the correct normalization is J ip 2 (r )d 3 r — 1 and hence ^ e (0) oc E z l 2 , 
where E is energy of the electron. There are two terms in the right hand side of eq. ( Al ) , the first term is proportional 
to E z l 2 and the second one is propo rtio nal to E 9 / 2 . For an external electron E — > and hence the second term must be 
neglected. Afte r that the equation (Al) is getting linear in t/j and we can return to the normalization ipg(0) — 3/(87r). 
Similar to (A5) it is convenient to use the substitution 



where f e satisfies the equation 



He(r) = \j 7^ Z (P ■ PN)fe(p), 



A„ + - / e (p) = -4tt%). 



(A10) 



(All) 



Note that due to the different normalization of ip e the coefficient in the right hand side of eq. (A10), including power 
of Z is different from that in eq. (A5). Solution of eq. (All) is 
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(A12) 



where Ni(x) is the singular Bessel function. At small x this function behaves as Ni ~ —2/irx, therefore f e (p) has 
correct behavior at small p: f e (p) ~ 1/p- Together with ( A10) this gives 



OiPe = ~\ -Z-Z{H ■ fJ, N ) . 

V 3 sfp 



(A13) 



To make sure that this is the correct solution one has to prove validity of the orthogonality condition ( A4). With eqs. 
(§) and (|A13|) one finds that the overlapping is of the form 

(A14) 



<<ty„|^„>oc / J 1 W8p)N 1 (^8p)pdp k / J^N^x'dx. 
Jo Jo 

This integral is not well defined at oo. The origin for this is clear: we are working with zero energy state. To correct 

rt 1 

the situation, one has to introduce an exponential factor e~ px and then to consider the limit p — > 0. 



/>00 -i -i 

(<ty«|Vn) « lim / e-^Mx^ix^dx cx lim — e -^W 3 72,3/2(1//?) « Mm -^e^K 3 



p^o 



0^0 P 5 



p^o f3 13 / 2 



0. 



(A15) 



For the evaluation of the integral we have used ref. |16y. Formula (A15) completes the prove of the orthogonality. 
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